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ABSTRACT 

Recent in situ Ulysses and Galileo observations of the source regions of type III solar radio bursts appear to 
show an absence of ion acoustic waves 5 produced by nonlinear Langmuir wave processes such as the electrostatic 
(ES) decay, in contradiction with earlier ISEE 3 observations and analytic theory. This Letter resolves these 
apparent contradictions. Refined analyses of the maximum 5-wave electric fields produced by ES decay and of 
the characteristics of the Ulysses Wave Form Analyzer (WFA) instrument show that the bursty 5 waves observed 
by [SEE 3 should be essentially undetectable by the Ulysses WFA. It is also shown that the maximum 5-wave 
levels predicted for the Galileo event are S the instrumental noise level, thereby confirming an earlier suggestion. 

Thus, no contradictions exist between the ISEE 3 and Ulysses! Galileo observations, and no evidence exists against 
ES decay in the published Ulysses and Galileo data. All available data are consistent with, or at worst not 
inconsistent with, the ES decay proceeding and being the dominant nonlinear process in type III bursts. 

Subject headings: interplanetary medium — plasmas — solar wind — Sun: radio radiation — waves 
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1. INTRODUCTION 

Over 30 years ago Ginzburg & Zheleznyakov (1958) sug- 
gested that nonlinear Langmuir wave processes are important 
in type III solar radio bursts. The processes considered since 
then include scattering off thermal ions (e.g., Ginzburg & 
Zheleznyakov 1958; Melrose 1980), three-wave coalescence 
and decay processes (e.g., Melrose 1982; Lin et al. 1986; 
Robinson, Willes, & Cairns 1993a; Thejappa et al. 1993a; 
Robinson, Cairns, & Willes 1994; Cairns & Robinson 1995) 
and modulational instabilities (e.g., Papadopoulos, Goldstein, 
& Smith 1974; Thejappa et al. 1993a, b). Questions remain as 
to the identities and roles of the important Langmuir wave 
processes. For instance, although scattering off thermal ions is 
well known to be ineffective, it is unresolved which three-wave 
and modulational (four-wave) processes are active and 
whether these processes quantitatively limit the Langmuir 
wave growth, allow persistence of the electron beam to 
distances ~1 AU, and generate the observed radiation. In- 
deed, comparisons of the recent stochastic growth theory with 
ISEE 3 data (Robinson 1992; Robinson, Cairns, & Gurnett 
1992, 1993b; Robinson et al. 1993a, 1994) imply that (linear) 
stochastic growth effects alone are responsible for limiting the 
growth of all but the largest amplitude Langmuir waves, for 
the observed wave clumping, and for the beam’s persistence. 
These authors argue that specific three-wave processes are 
responsible for limiting the highest amplitude Langmuir waves 
and for generating the radiation. 

The existence (or not) and frequencies of low-frequency ion 
acoustic waves produced in association with intense Langmuir 
waves, or fine structure in the Langmuir spectra and wave- 
forms, constitute the simplest observational routes for testing 
which nonlinear processes are active in type III bursts. Lin et 
al. (1986) reported ISEE 3 observations near 1 AU of such 
product waves. They found the wave frequencies to be 
commensurate with those expected (wavenumbers ^Wp/t/j,) 
for the electrostatic (ES) decay L -^L' + S and the electro- 
magnetic (EM) decay L T{fp) + 5'. Here L and L' de- 


note primary and product Langmuir waves, respectively, 5 and 
S' are ion acoustic waves, T{fp) is a transverse EM wave with 
fi-equency close to the electron plasma frequency /p, Wp = litfp, 
and Vi, is the beam speed. Lin et al. (1986) argued against 
modulational instabilities being relevant and found the wave 
frequencies to be most consistent with the EM decay (despite 
analytic theory definitely favoring the ES decay). However, a 
refined analysis of these same data (Cairns & Robinson 1995) 
found strong evidence that the ES decay produces the great 
majority of these waves. The data are also consistent with the EM 
decay proceeding (and producing radiation near/p) when stim- 
ulated by the ES decay. In fact, three lines of evidence support 
the ES decay being the dominant nonlinear process in type III 
bursts: the frequencies and onset times of the observed low- 
frequency waves (Cairns & Robinson 1995), a cutoff in the 
probability distribution P(£l) at high L-wave electric fields El 
near the predicted threshold for ES decay (Robinson et al. 1993a, 
b), and the frequencies of beat waveforms observed by the 
Galileo spacecraft in a type III source region (Gurnett et al. 
1993). The first two lines of evidence are derived from ISEE 3 
data. 

At first sight, Ulysses observations from 1 to 4 AU appear 
totally at variance with this picture. Thejappa et al. (1993a) 
discussed three type III bursts with local radio emission and L 
waves but with no 5-wave products above instrumental back- 
ground. They calculated that the maximum 5-wave fields 
predicted for ES decay should be well above the noise level 
and similar to the ISEE 3 observations and thus argued that 
the absence of decay products implies that the ES and EM 
decays do not proceed, thereby contradicting the ISEE 3 
observations. Thejappa et al. (1993a, b) then argued that 
modulational instabilities produce, and the dominant radia- 
tion processes only involve, the localized Langmuir wave 
packets observed. However, no direct evidence for the modu- 
lational instability exists due to a lack of product waves. 

Galileo observations of a type III burst near 1 AU (Gurnett 
et al. 1993) showed Langmuir waveforms with beats at fre- 
quencies consistent with ES decay, as predicted by Cairns & 
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Robinson (1992), but found no direct evidence for 5-wave 
products. Gumett et al. (1993) suggested that this might 
simply be a result of high noise levels but did not analyze the 
situation further. 

This Letter investigates and resolves the apparent contra- 
dictions between the ISEE 3 and UlyssesIGalileo data, address- 
ing the very different interpretations adopted for the ISEE 3 
and Ulysses observations. Sections 2 and 3 present refined 
predictions for the maximum 5-wave fields produced by ES 
decay for the Ulysses events, taking into account the charac- 
teristics of the relevant Ulysses wave instrument and using in 
situ plasma data. It is shown that 5-wave bursts like those 
observed by ISEE 3 should be essentially undetectable in the 
time-averaged Ulysses data. Analogous calculations in § 4 for 
the Gurnett et al. (1993) Galileo observations show that the 
maximum predicted 5-wave fields for ES decay are less than or 
similar to the typical noise level. The key conclusion is that the 
apparent contradictions between the ISEE 3 data and the 
Ulysses and Galileo data do not exist. 


2. THEORETICAL FRAMEWORK 

The ES decay L-^L' -F 5 proceeds fastest as a process that 
scatters a beam-driven L wave, with ~ utpvjvl, into a 
backward going L' wave (with ~ — kz.) and a forward-going 

5 wave with ks 2ki_ (e.g., Melrose 1982; Cairns 1987). 
Frequency and wavevector conservation and wave dispersion 
relations constrain the wave frequencies and wavevectors. A 
standard relation exists between the electric field energy 
density and the occupation number Ny of quanta in mode M. 


io 
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( 1 ) 


where Ruiku) is the ratio of electric to total energy in the 
mode M. Considering conservation of the total number of 
quanta interacting and the volumes of wavevector phase space 
occupied by the participating waves, the maximum 5-wave 
field is related semiquantitatively to the L-wave field by 
(Robinson et al. 1993a, eq. [47c]) 


Es 



( 2 ) 


Here Vg is the ion acoustic speed defined by F| = ksT^{\ + 3TJ 
T^)lnii, and and T, are the electron and ion temperatures. 

The limit of equation (2) is lower by a factor of 2 than 
the Thejappa et al. (1993a, p. 842) prediction. The factor of 2 
refinement is caused by geometry, the solid angles Aflj and 
Aflz, occupied by the 5 waves and primary L waves are related 
by Aflg « Alii, (Robinson et al. 1993a) and not by the 
TThejappa et al. relation Aflj = 4Afl/,. The factor of 2^'^ 
change arises from requiring conservation of the total number 
of quanta (Robinson et al. 1993a), rather than assuming 
conservation of occupation numbers. 


3. t/Z,ySS£S ANALYSIS 

Thejappa et al. (1993a) analyzed Ulysses observations of 
three type III bursts which exhibited local radio emission and 
Langmuir waves: one each in the periods 1300-1500 UT on 
1990 December 11, 0630-0800 UT on 1991 February 22, and 
1700-2000 UT on 1991 March 7. In each case S 10 kHz. 
Two Ulysses detectors are relevant here (Stone et al. 1992): the 
Fast Envelope Sampler (FES) and the Wave Form Analyzer 
(WFA). For these bursts, the FES provided the L-wave 


TABLE 1 

Ulysses Measurements of and T, 


Range of T, Range of T, T, T, 

Event (K) (K) (K) (K) 


Dee II 7.3-11 x 10“ 2.S-6.4 X 10“ 1.0 X 10^ 4.5 x 10“ 

Feb 22 7.4-7.9 X 10“ 1.0-1.7 X 10“ 7.6 X 10“ 1.3 x 10“ 

Mar 7 4.0-4.9 X 10“ 0.7-1.9 x 10“ 4.5 X 10“ 1.3 X 10“ 


envelope fields every 1 ms in blocks 1024 samples long in a 
6-60 kHz bandwidth (the detector time constant is less than 1 
ms). Owing to telemetry restrictions, only the FES period with 
the highest peak field is retained and reported approximately 
every half-hour. Detailed operation of the WFA is described 
below; for now, we note that the WFA covers the frequency 
range for ES decay product waves and provides outputs once 
every DPU cycle (64 or 128 s, depending on the telemetry rate) 
with a large amount of instrumental dead time. Since both the 
L and 5 waves observed by ISEE 3 (Lin et al. 1986) and the L 
waves observed by Ulysses (Kellogg et al. 1992; Thejappa et al. 
1993a) are extremely spiky, the slow sampling rate and large 
dead time of the \\EA instrument makes it very difficult to 
detect ES decay products, as we show below. 

Two further refinements to the previous analysis are to use 
in situ Ulysses measurements for T, and T, in equation (2) 
rather than nominal values, and to consider time variations in 
vt, (e.g., Lin et al. 1981, 1986) rather than a single nominal 
value. In fact, the measured and T, are substantially lower 
than the Thejappa et al. (1993a) assumed values T, = 2 X 10^ 
K and T, = 4 X 10^ K. The effect of this underestimation can 
be seen by rearranging equation (2) as 

Es s ' (1 + 3r,./rj^«, (3) 

where F, = is the electron thermal speed. Reduc- 

ing T, or the ratio T,/r, therefore reduces the maximum 
5-wave field Eg^ predicted for ES decay products. 

Table 1 summarizes Ulysses data for and T-, from the 
SWOOPS instrument for the three events (J. L. Phillips 1995, 
private communication; see Phillips et al. 1993), giving ob- 
served ranges and typical values of T, and 7). The quoted T, 
values are those for the “core” component of the solar wind 
electron distribution, since Fs is small compared with the core 
and halo thermal speeds. The quoted T, are averages of the 
maximum and minimum values measured. Since Ulysses does 
not measure the beam velocity, we consider a range u* = 3.5- 
10 X 10’ m s“' based on ISEE 3 data (Lin et al. 1981, 1986; 
Cairns & Robinson 1995). 

The maximum ratios of 5-wave to L-wave fields calculated 
from equation (3), based on Table I’s typical T, values and 
the range of vg, are given in Table 2’s third column. These 
ratios are then multiplied by the maximum FES L-wave fields 
(col. [4]) to give Limax- The quoted FES fields correspond to 
the broad wave packets seen, not the superposed 1 ms spikes 
sometimes seen; the 1 ms spikes have spatial scales ~10 times 
too small to excite ES decay products in the range 100-300 Hz 
range observed by ISEE 3. Column (5) gives the value of Lj„,ax 
observable by an ideal FES-like detector. As described more 
below, instrumental limitations and the intrinsic burstiness of 
the L and 5 waves limit the detectability of these 5-wave fields. 
The final column contains the noise level in the WFA 112 Hz 
“channel” for each event, obtained from the Thejappa et al. 
(1993a) figures as the square root of the background spectral 
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TABLE 2 

Maximum Predicted 5-Wave Levels and the Ulysses Instrumental Noise Levels 


vb Ratio FES £(, Ideal Es Noise 

Event (iHs"’) EsIEl (mV m“') (jiV m“') (/xV m“') 

0) (2) (3) (4) (5) (6) 

Dec 11 10 X 10’ 1.1 X 10“’ 3.0 3.4 2.8 

3.5 X 10’ 5.4 X 10"’ 3.0 16 2.8 

Feb 22 10X10’ 6.5 x lO-"* 2.5 1.6 1.0 

3.5 X 10’ 3.1 X 10"’ 2.5 7.9 1.0 

Mar 7 10 X 10’ 5.2 X 10“' 5.0 2.6 1.2 

3.5 X 10’ 2.5 X 10“’ 5.0 12 1.2 


density times the channel bandwidth. Table 2 shows that 
^srnax S 16 ixV m“', less than 2-10 times the noise background. 
These fields are smaller than the Thejappa et al. (1993a) 
quoted upper limit of 32 /xV m“’ owing to the use of measured 
plasma parameters and the refined prediction in equation (3). 

Before detailed consideration of the WFA’s characteristics 
(Stone et al. 1992; P. J. Kellogg & R. J. MacDowall 1995, 
private communications), it might appear that ES decay 
products should be easily detectable. This is, however, not the 
case. The WFA’s input is a bandpass-limited signal sampled at 
1792 FIz. The WFA measures a spectrum for each octant of 
spacecraft spin (spin period ~12 s), with spectral densities 
calculated at specific frequencies (448, 224, 149, 112 Hz, etc.) 
using windowed digital fourier transforms on four cycle 
lengths (16 samples) of the field time series (appropriately 
time-averaged). The peak spectrum for each octant is con- 
structed using the maximum spectral density at each “channel” 
frequency every DPU cycle (64 or 128 s). Average spectra are 
constructed for each “quadrant” of spacecraft spin by averag- 
ing over opposing octants and then over the spin periods in 
each DPU cycle. For these bursts, the WFA reported four 
time-averaged quadrant spectra and eight peak octant spectra 
every 64 s. The analysis proceeds in four stages: 

1. Using only a four-cycle time interval every 1.5 s for 
calculating the spectral density at each frequency means 
that the instrument has substantial dead time at high 
frequencies (P. J. Kellogg 1995, private communication). 
For instance, the dead time at 112 Hz is (1.5-4/112) s, or 
~98% of the time between spectra; the percentage of dead 
time decreases at the lower frequencies. The probability 

of the WFA sampling at a given time is thus ~2% at 112 Hz 
or above, as a result of the long instrumental dead time. 

2. Analysis of the Cairns & Robinson (1995) detailed ISEE 3 
L- and 5-wave data indicates that the probability, Pbursi. of 

f an 5-wave burst produced by ES decay occurring at a given 

time ranges from ~1% to 3% at 1 AU. 

J 3. The probability of the WFA being on and sampling an 
5-wave spike is then the product PonPburst ~0.02%-0.06% at 
1 AU. Therefore the WFA needs to process 2000-5000 
spectra on average in order to sample one 5-wave burst. 
Since the WFA calculates 40 octant spectra every 64 s, this 
number of spectra corresponds to an event duration of 
50-130 minutes. That is, the WFA is predicted to sample 
only about one 5-wave burst every 1-2 hr; this corresponds 
to approximately one 5-wave burst for each of Thejappa et 
al.’s (1993a) events. 

4. Thejappa et al. (1993a) present WFA data from the time- 
and quadrant-averaged spectra. Consider the 40 octant 
spectra used to construct this 64 s average and suppose that 


at some frequency n spectra are enhanced by an average 
factor a above the true noise level A,. Then the average is 


n(a — 1) 


Table 2’s last two columns imply that a S 10 for ES decay 
products. For n = 1 (justified by point 3 above), then 
xljpp « 1.L4, for a = 5, and A^pp « 1.2x4, for a = 10. Thus, 
even if an ES decay spike occurs in one 64 s period it will 
increase the apparent noise level by S20%. Enhancements 
of this magnitude are essentially undetectable in Thejappa 
et al.’s (1993a) data owing to the presence of multiple, fairly 
periodic variations in the spacecraft’s average noise level of 
similar or larger size. 

In summary, the mode of operation of the Ulysses WFA 
instrument and the rarity and burstiness of the ES decay 
products observed by ISEE 3 lead to the prediction that 
5-wave signals should be extremely rare in the Ulysses WFA 
data and essentially unobservable in the average WFA output. 
Thejappa et al.’s (1993a) result of an apparent absence of 
5-wave signals in the WFA data is therefore theoretically 
expected. The best chance of observing ES decay products 
with Ulysses is in the peak WFA spectra. Even in the peak 
data, however, the maximum 5-wave fields predicted are only 
~2-10 times the average noise level, with one 5-wave detec- 
tion expected every 1-2 hr. An extremely low likelihood exists 
then of a successful detection, consistent with Thejappa et al. 
(1993a) not identifying 5-wave events in the peak data. The 
published Ulysses data therefore provide no evidence for or 
against the ES decay proceeding, there is no contradiction 
between the Ulysses and ISEE 3 observations, and there is also 
no necessity (based on these data) to appeal to modulational 
instabilities to limit L-wave growth and reduce the beam’s 
energy losses. The present Ulysses data are therefore not 
inconsistent with theories, based on ISEE 3 data that rely 
heavily on ES decay (e.g., Cairns & Robinson 1992, 1995; 
Robinson et al. 1993a, b, 1994). 


4. Gx4L/L£0 ANALYSIS 

Gurnett et al. (1993) reported Galileo observations at 0.98 
AU of the L waves, radio emission, and electrons of a type III 
burst on 1990 December 10. Beats in the Langmuir waveforms 
had frequencies consistent with ES decay, as predicted by 
Cairns & Robinson (1992). However, no 5-wave products of 
ES decay were observed directly; Gurnett et al. (1993) sug- 
gested this was the result of high instrumental noise levels. 
This suggestion can be evaluated similarly to the Ulysses 
observations. Gurnett et al. (1993) presented wideband re- 
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TABLE 3 


Maximum Predicted S-Wave Levels and the Galileo Instrumental Noise Levels 


Vb 

Ratio 

Wide-band Ei_ 

Ideal Es 

Channel 

Max Es 

100 Hz 

316 Hz 

(ms ') 

EsIEl 

(mV m“') 

(pV m-') 

(mV m“‘) 

(jiV m“‘) 

(ftV m-') 

(nV m“‘) 

(1) 

(2) 

( 3 ) 

(4) 

( 5 ) 

( 6 ) 

( 7 ) 

(8) 

10 X loP .... 

. 1.3 X 10-’ 

1.7 

2.2 

0.5 

0.6 

4.5 

2.3 

3.5 X 10’ ... 

. 6.1 X 10-’ 

1.7 

10 

0.5 

3.1 

4.5 

2.3 


ceiver and spectrum (channel) analyzer data for the L waves, 
but only the channel analyzer sampled the frequency range 
predicted for ES decay products. 

Two important effects influence predictions for the maxi- 
mum S-wave fields observable by Galileo. First, comparison of 
Gumett et al.’s (1993) wide-band and channel L-wave data 
(their Figs. 6 and 2, respectively) implies that the 50 ms 
detector time constant and slow sampling rate (once every 37.3 
s) for the channel fields cause them to be underestimated by a 
factor ~3.4 (= 1.7 mV m“V0.5 mV m"’) compared with the 
maximum fields seen by the wide-band receiver. The S-wave 
fields are sampled using the low-frequency channels, also with 
time constants ~50 ms, thereby being subject to the same 
underestimation when compared with predictions based on 
the wide-band L-wave fields. The maximum observable 5-wave 
field predicted, Esm^^, is therefore reduced by the same 
factor of 3.4 compared with a prediction that does not take 
into account time constant and sampling rate effects. 

Second, Gumett et al. (1993) used a nominal T, = 1.3 X 10^ 
K as a result of the Galileo thermal plasma instmment being 
switched off. Taking other nominal parameters TJTi = 3 and 
Vb = 3-10 X 10’ m s~‘ and including the underestimation 
factor, the predictions for ^imaxoiw are given in Table 3’s 
column (6): E^maxobs ~ 1-3 m^V m"‘. (Col. [4] shows the ideal 
Esmsx, i e-> without the underestimation factor.) 

The instmmental background fields were investigated using 
averages of four spectrum analyzer channels (with intrinsic 
bandwidths ~5 Hz) centered at 100 Hz and 316 Hz, with time 
resolution of 32 s provided by D. A. Gumett & L. J. Granroth 
(1995, private communications). For the period 0800-0830 
UT, before the L-waves started, the typical 100 Hz fields were 
~4.5 /xV m”' with minimum and maximum fields of 2.8 and 24 
/xV m“‘. The typical 100 Hz fields remained at levels ~5 /xV 
m'’ for the event period 0830-1000 UT, while the minimum 
and maximum fields were 1.5 and 65 /xV m"‘. At 316 Hz the 
fields showed little variation over the entire period 0800-1000 
UT: the typical field remained ~2.3 /xV m“‘, and the mini- 
mum and maximum fields observed were 1.5 and 4.0 fxV m"'. 
Table 3 makes it clear that Ei„ax.obs is predicted to be below or 
similar to the typical noise level and only slightly above the 
lowest channel fields observed. We therefore confirm the 
Gumett et al. (1993) suggestion that sensitivity effects 
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should prevent reliable detection of ES decay products for 
this event. (Also, Galileo’s low sampling rate implies that it 
should sample S-waves ~10 times less often as ISEE 3.) 

5. CONCLUSIONS 

The detectability of ion acoustic products of the ES decay 
L^L' + S has been examined in detail for recent Ulysses and 
Galileo observations of type III burst source regions. In both 
cases, the nonobservation of ES decay products has been 
explained: waves similar to those observed by ISEE 3 (Lin et 
al. 1986; Cairns & Robinson 1995) could be present but would 
be essentially undetectable using the Ulysses or Galileo plasma 
wave detectors. For the Ulysses events discussed by Thejappa 
et al. (1993a) this nondetectability is due to (1) the maximum 
5-wave fields being less than a factor of 10 above the back- 
ground, (2) the use of quadrant; and time-averaged WFA 
data, which should limit the effect of the maximum 5-wave 
burst to raising the nominal noise level by S20% (thereby 
rendering bursts essentially undetectable), and (3) the large 
instmmental dead time of the WFA and the burstiness of the 
(ISEE 3) 5 waves, implying that the WFA should very rarely 
detect an 5-wave burst (only one every 1-2 hr on average). The 
lack of ES decay products in the Galileo data is due to the 
predicted maximum 5-wave fields being less than or of order 
the typical noise level, thereby confirming Gumett et al.’s 
(1993) suggestion. The available Ulysses WFA and Galileo 
channel data thus provide no evidence for or against the ES 
decay proceeding, there is no contradiction between the ISEE 
3 observations and the Ulysses and Galileo observations, and 
these data do not argue that modulational instabilities limit 
the Langmuir wave growth (in the apparent absence of ES 
decay) for the Ulysses and Galileo events. Indeed, now all the 
available data are consistent with, or at worst not inconsistent 
with, the ES decay proceeding and being the dominant non- 
linear process in type III bursts, as predicted by analytic theory. 
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